Thermodynamics of enhanced heat transfer: a model study 
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Situations where a spontaneous process of energy or matter transfer is enhanced by an external 
device are widespread in nature (human sweating system, enzyme catalysis, facilitated diffusion 
across bio-membranes, industrial heat exchangers). The thermodynamics of such processes remains 
however open. Here we study enhanced heat transfer by a model junction immersed between two 
thermal baths at different temperatures Th and Tc {Th > Tc). The transferred heat power is 
enhanced via controlling the junction by means of external time-dependent fields. Provided that 
the spontaneous heat flow process is optimized over the junction Hamiltonian, any enhancement 
of this spontaneous process does demand consumption and subsequent dissipation of work. The 
efficiency of enhancement is defined via the increment in the heat power divided over the amount 
of consumed work. We show that this efficiency is bounded from above by Tc/{Th — Tc). Formally 
this is identical to the Carnot bound for the efficiency of ordinary refrigerators which transfer heat 
from cold to hot. It also shares some (but not all) physical features of the Carnot bound. 

PACS numbers: 05.70.-a, 05.30.-d, 05.70.Ln, 84.60.-h 



I. INTRODUCTION 

Heat spontaneously flows from hot to cold. In some 
cases it is necessary to reverse this flow. Devices that per- 
form this operation need an external input of high-graded 
energy (work), which is lost in the process: refrigerators 
cool a colder body in the presence of a hotter environ- 
ment, while heaters heat up a hot body in the presence 
of a colder one [1]. The efficiency (or coefficient of perfor- 
mance) of these devices is naturally defined as the useful 
effect — for refrigerators this is the heat extracted from 
the colder body, while for heaters this is the heat deliv- 
ered to the hotter body — divided over the work consumed 
per cycle from the work-source [1]. The first and second 
laws of thermodynamics limit this efficiency from above 
by the Carnot value: For a refrigerator (heater) operat- 
ing between two thermal baths at temperatures Tc and 
Th, respectively, the Carnot efficiency reads [1] 
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There are however situations, where the spontaneous 
direction of the process is the desired one, but its power 
has to be increased. An example of such a process is 
perspiration (sweating) of mammals [2] . A warm mam- 
malian body placed in a colder environment will nat- 
urally cool due to spontaneous heat transfer from the 
body surface. Three spontaneous processes are involved 
in this: infrared radiation, conduction and convection 
[2]. When the environmental temperature is not very 
much lower than the body temperature, the spontaneous 
processes are not sufficiently powerful, and the sweating 
mechanism is switched on: sweating glands produce wa- 
ter, which during evaporation absorbs latent heat from 
the body surface and thus cools it [2]. Some amount 
of free energy (work) is spent in sweating glands to wet 
the body surface. Similar examples of heat transfer are 
found in the field of industrial heat-exchangers, where 



the external source of work is employed for mixing up 
the heat- exchanging fluids. 

The main feature of these examples is that they com- 
bine spontaneous and driven processes. Both are macro- 
scopic, and with both of them the work invested in en- 
hancing the process is ultimately consumed and dissi- 
pated. Pertinent examples of enhanced transport exist in 
biology [4, 5]. During enzyme catalysis, the spontaneous 
rate of a chemical reaction is increased due to interac- 
tion of the corresponding enzyme with the reaction sub- 
strate. (A chemical reaction can be regarded as particle 
transfer from a higher chemical potential to a lower one.) 
There are situations where enzyme catalysis is fueled by 
external sources of free energy supplied by co-enzymes 
[4]. However, many enzymes function autonomously and 
cyclically: The enzyme gathers free energy from binding 
to the substrate, stores this free energy in slowly relaxing 
conformational degrees of freedom [6, 7], and then em- 
ploys it for lowering the activation barrier of the reacion 
thereby increasing its rate [4-7]. Overally, no free en- 
ergy (work) is consumed for enhancing the process within 
this scenario. Similar situations are realized in transport- 
ing hydrophilic substances across the cell membrane [4]. 
Since these substances are not soluble in the membrane, 
their motion along the (electro-chemical) potential gradi- 
ent is slow, and special transport proteins are employed 
to enhance it [4, 5]. Such a facilitated diffusion normally 
does not consume free energy (work). 

These examples of enhanced processes motivate us to 
ask several questions. Why is that some processes of en- 
hancement employ work consumption, while others do 
not? When enhancement does (not) require work con- 
sumption and dissipation? If the work-consumption does 
take place, how to define the efficiency of enhancement, 
and are there bounds for this efficiency comparable to 
(1)? These questions belong to thermodynamics of en- 
hanced processes, and they are currently open. Laws 
of thermodynamics do not answer to them directly, be- 



cause here the issue is in increasing the rate of a process. 
Deahng with time-scales is a weak-point of the general 
thermodynamic reasoning [3], a fact that motivated the 
development of finite-time thermodynamics [8] . 

Here we address these questions via analysing a quan- 
tum model for enhanced heat transfer. The model de- 
scribes a few-level junction immersed between two ther- 
mal baths at different temperatures; see section II. The 
junction is subjected to an external field, which enhances 
the heat transferred by the junction along its spontaneous 
direction. The virtue of this model is that the optimiza- 
tion of the transferred heat over the junction Hamiltonian 
can be carried out explicitly. Based on this, we determine 
under which conditions the enhancement of heat-transfer 
does require work-consumption. We also obtain an upper 
bound on the efficiency of enhancement, which in several 
aspects is similar to the Carnot bound (1). 

Heat flow in microscale and nano-scale junctions re- 
ceived much attention recently [9-16, 19]. This is re- 
lated to the general trend of technologies towards smaller 
scales. Needless to stress that thermodynamics of en- 
hanced heat-transfer is relevant for this field, because it 
should ultimately draw the boundary beetween what is 
possible and what is not when cooling a hot body in the 
presence of a colder one. Brownian pumps is yet another 
field, where external fields are used to drive transport; 
see, e.g., [20, 21] and references therein. Some of the set- 
ups studied in this field are not far from the enhanced 
heat transport investigated here. However, thermody- 
namical quantities (such as work and enhancement effi- 
ciency) were so far not studied for these systems, though 
thermodynamics of Brownian motors [work-extracting 
devices] is a developed subject reviewed in [22]. 

The rest of this paper is organized as follows. The 
model of heat-conducting junction is introduced in sec- 
tion II. Section HI shows how the transferred heat (with 
and without enhancing) can be optimized over the junc- 
tion structure. The efficiency of enhancing is studied in 
section IV. Section V discusses how some of the obtained 
results can be recovered from the formalism of linear non- 
equilibrium thermodynamics. We summarize in section 
VI. Several questions are relegated to Appendices. 



II. THE MODEL 

Our model for the heat pump (junction) consists of 
two quantum systems H and C with Hamiltonians Hn 
and He, respectively; see Fig. 1. Each system has n 
energy levels and couples to its thermal bath. Similar 
models were employed for studying heat engines [17, 18] 
and refrigerators [19]. It will be seen below that this 
model admits a classical interpretation, because all the 
involved initial and final density matrices will be diago- 
nal in the energy representation. We shall however work 
within the quantum framework, since it is more intuitive. 

Initially, H and C do not interact with each other. Due 
to coupling with their baths they are in equilibrium at 
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FIG. 1: The heat pump model. The few-level systems H 
and C operate between two baths at temperatures Tc and Th 
Tc < Th- During the first step of operation the two systems 
interact together either spontaneously or driven by a work- 
source at the cost of work W. During this stage couplings 
with the thermal baths is neglected (thermal isolation). In 
the second step the systems H and C do not interact with 
each other and freely relaxes to their equilibrium states (2) 
under action of the corresponding thermal bath. 

temperatures T/j — 1/ (ih > Tc = I//3c [we set ks ~ !]■ 
^ = g-/3;.ffH/tj.[e-''''^«], a = e-'^''"^/tT[e-^^"^], (2) 

where p and a are the initial Gibbsian density matrices 
of H and C, respectively. We write 



P = diag[' 



■,ri\ 



CT = diag[s„,...,si], (3) 



Hh = diag[e„, ...,ei = 0], He = diag[^„, ...,^i = 0], 

where diag[a,..,&] is a diagonal matrix with entries 
(a,..., 6), and where without loss of generality we have 
nullified the lowest energy level of both H and C . Thus 
the overall initial density matrix is 

f^in ^ P®(T, (4) 

and the initial Hamiltonian of the junction is 

i?o = i?H®l + l®i?c. (5) 

A. Spontaneous operation 

During a spontaneous process no work is exchanged 
with external sources. For our situation a spontaneous 
heat transfer will amount to a certain interaction be- 
tween H and C. Following to the approach of [24-26] 
we model this interaction via a Hamiltonian that con- 
serves the (free) Hamiltonian Hq [see (5)] for all interac- 
tion times. This then realizes the main premise of spon- 
taneous processes: no work exchange at any time. Our 
model for spontaneous heat transfer consists of two steps. 

1. During the first step H and C interact with each 
other [collision]. We assume that this interaction takes 
a sufficiently short time S, and during this time the cou- 
pling with the two thermal baths can be neglected [ther- 
mal isolation] . The interaction is described by the Hamil- 
tonian Hint added to (5): 
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The overall Hamiltonian H again lives in the n^- 
dimensional Hilbcrt space of the junction ^. 

As argued above, the interaction Hamiltonian com- 
mutes with the total Hamiltonian: 



[Ho, Hint] — 0, 



(7) 



making the energy Hq a conserved quantity ^ . To have a 
non-trivial effect on the considered system, the interac- 
tion Hamiltonian Hmt should not commute with the sep- 
arate Hamiltonian: [Hh (g) 1, iJjnt] 7^ 0. For this to be the 
case the spectrum of Hq should contain at least one de- 
generate eigenvalue. Otherwise, relations [ifo,-ffint] — 
and [Hu (g) 1, ifo] =0 will imply [Hu ® 1, Hint] = (and 
thus a trivial effect of Hmt), because the eigen-base of Hq 
will be unique (up to re-numbering of its elements and 
their multiplication by phase factors). The energy 



Q^l = tr ( Fh 



p- trc 



Irt IT 

-,--i--nii 



fl.n ei 



Hi, 



, (8) 



lost by H during the interaction is gained by C. Here 
trn and trc sue the partial traces. 

Commutative interaction Hamiltonians (7) are applied 
to studying heat transfer in [24-26]. Refs. [24, 25] are 
devoted to supporting the thermodynamic knowledge via 
quantum Hamiltonian models. In contrast, the approach 
of [26] produced new results. 

2. For times larger than S, H and C do not interact 
and freely relax back to their equilibrium states (2, 4) 
due to interaction with the corresponding thermal baths. 
These equilibrium states are reached after some relax- 
ation time T. Thus the cycle is closed — the junction re- 
turns to its initial state — and Qf^ given by (8) is the 
heat per cycle taken from the hot thermal bath during 
the relaxation (and thus during the overall cycle). 

It should be obvious that once T^ > Tc we get Qj^ > 
0: heat spontaneously flow from hot to cold. The proof 
of this fact is given in [18, 19, 24-26]. 

For times larger than r there comes another interaction 
pulse between H and C, and the cycle is repeated. 



1. Power 

Recall that the power of heat-transfer is defined as 
the ratio of the transferred heat to the cycle duration r. 



^ More precisely, we had to write the Hamiltonian (6) as Hu 8i 1 + 
1 ® He + Q(t)ii'int , where a{t) is a switching function that turns 
to zero both at the initial and final time. This will however not 
alter the subsequent discussion in any serious way. 

^ This implementation of spontaneous heat-transfer processes ad- 
mits an obvious generalization: one need not require the conser- 
vation of Ho for all interaction times, it suffices that no work 
is consumed or released within the overall energy budget of the 
process in the time-interval [0, S]. For our purposes this general- 
ization will not be essential; see (27). 



Q^ /t. For the present model r is mainly the duration 
of the second stage, i.e., r is the relaxation time, which 
depends on the concrete physics of the system-bath cou- 
pling. For a weak system-bath coupling t is larger than 
the internal characteristic time of H and C. In contrast, 
for the collisional system-bath interaction, r can be very 
short; see Appendix A. Thus the cycle time t is finite, 
and the power Ofj^ /t does not vanish due to a large 

cycle time, though it can vanish due to Qfi^ — > 0. 

Note that some entropy is produced during the free 
relaxation. This entropy production can be expressed 
via quantities introduced in (4-8); see [19] for details. 
The first step does not produce entropy, because it is 
thermally isolated and is realized by a unitary operation 
that can be reversed by operating only on observable de- 
grees of freedom (H 4- C). It is seen that no essential 
aspect of the considered model depends on details of the 
system-bath interaction. This is an advantage. 



B. Driven operation 

The purpose of driving the junction with an external 
field is to enhance [increase] the spontaneous heat Q}^ ■ 
The driven regime reduces to the above two steps, but 
instead of the spontaneous interaction we have the fol- 
lowing situation: the interaction between H and C is 
induced by an external work-source. Thus (7) does not 
hold anymore. The overall interaction [between H, C 
and the work-source] is described via a time-dependent 
potential V{t) in the total Hamiltonian 



Hii®l + l®Hc + V{t) 



(9) 



of H + C. The interaction process is still thermally 
isolated: V{t) is non-zero only in a short time- window 
< t < (5 and is so large there that the infiuence of the 
couplings to the baths can be neglected. 

Thus the dynamics of H + C is unitary for < t < (5: 



f7f 



^(5)=U^iU\ U^Te-T^io^'^v{s)+H,]^ (^Q^ 



where Q,i — il(0) = p®a is the initial state defined in (2), 
ilf is the final density matrix, U is the unitary evolution 
operator, and where 7" is the time-ordering operator. The 
work put into H -(- C reads [1, 3] 

W ^Ef^Ei=ti[{H^®l + l® He) (^f - ^i) ], (11) 

where Ef and Ei are initial and final energies of H + C. 
Due to the interaction, H (C) looses (gains) an amount 
of energy Qh {Qc) 



Qh = tr{HH[p~ trc^f]), 
Qc^tr{Hc[tTH^f-cr])- 

Eqs. (11, 12) imply an obvious relation 



(12) 
(13) 



(14) 



Recall that for spontaneous processes not only the con- 
sumed work is zero, W ~ 0, but also the stronger con- 
servation condition (7) holds. 

Once H + C arrives at the final state ilun, the inter- 
system interaction is switched off, and H and C sep- 
arately [and freely] relax to the equilibrium states (2). 
During this process Qh is taken as heat from the hot 
bath, while Qc is given to the cold bath. Note from sec- 
tion II A 1 that the driven operation does not change the 
cycle time r, because the latter basically coincides with 
the relaxation time. Therefore, increasing Qh means in- 
creasing heat transfer power. 



III. MAXIMIZATION OF HEAT 
A. Maximization variables 

The type of questions asked by thermodynamics con- 
cerns limiting, optimal characteristics. Sometimes the 
answers are uncovered directly via the basic laws of ther- 
modynamics, an example being the Carnot bound (1). 
However, more frequently than not, this goal demands 
explicit optimization procedures [8] . 

We shall maximize the heat Qh transferred from the 
hot bath over the full Hamiltonian of the junction. For 
spontaneous processes this amounts to maximizing over 
Hamiltonian (6) living in the n^-dimensional Hilbert 
space of the junction H -I- C and satisfying condition (7). 
For driven processes we shall maximize over Hamiltoni- 
ans (9). In this case we shall impose an additional condi- 
tion that the work put into H + C in the first step does 
not exceed E > 0: 



W < E. 



(15) 



Once the work put into the system is a resource, it is 
natural to operate with resources fixed from above. 

Recall that the Hamiltonians (6, 9) live in the n^- 
dimensional Hilbert space. The bath temperatures Tc 
and Th and the dimension n^ (the number of energy lev- 
els) will be held fixed during the maximization. 

Due to (8) the maximization of the spontaneous heat 
Qh over the Hamiltonians (6, 7) amounts to maximiz- 
ing over unitary operators e^^''='*, and over the energies 
{^fc}fe=2' {^^k}k=2 °f H and C. Likewise, as seen from 
(9-11), the maximization of the driven heat Qh amounts 
to maximizing under condition (15) over all unitary op- 
erators U living in the n^-dimcnsional Hilbert space, and 
over the energies {£fe}^^2' {f^k}k=2- 

We should stress that the class of Hamiltonians liv- 
ing in the n^-dimensional Hilbert space [with or with- 
out constraint (7)] is well-defined due to separating the 
heat transfer into two steps (thermally isolated interac- 
tion and isothermal relaxation). More general classes of 
processes can be envisaged. For instance, we may write 
the free Hamiltonian as Hq + Hb,c + Hg^h, where Hq, 
Hbc and Hb.h are, respectively, the Hamiltonians of the 



junction and the two thermal baths. Now the Hamilto- 
nian iJint of spontaneous processes will be conditioned 
as [iJint, Hq + Hb.c + Hb.h] = 0. This condition is more 
general than (7) . Then the corresponding class of driven 
processes can be naturally defined via the same class of 
Hamiltonians but without this commutation condition. 
We do not consider here such general processes, since we 
are not able to optimize them. 



B. Unconstrainted maximization 

As seen below, the maximization of the spontaneously 
transferred heat (8) amounts to a particular case of max- 
imizing Qh- So we shall directly proceed to maximizing 
the driven heat Qh; see (12). 

First, take in (15) the simplest case: E ~ +oo. This 
case contains the pattern of the general solution. Here 
we have no constraint on maximization of Qh and it is 
done as follows. Since tr[Hiip] depends only on {sk}k=2^ 
we choose {lJ,k}k=2 ^^"^ ^(^) ^^ ^'^^^ ^he final energy 
tr[-ffH trci^f] attains its minimal value zero. Then we 
shall maximize tilHup] over {£k}k=2- Note from (3) 



£ji , . . . , £ji J , 

• 1 ^n^l 1 ■ ■ • 7 ^n^n I 



Hu<S)l = diag[ei ,..., ei,.. 
fli = p(g)a = diag[risi,...,ris„. 

It is clear that tr [iJn trcf^f] = ti [{Hh ^ l)ll ^iH'' ] 
goes to zero when, e.g., S2 = . . . = s„ ^^ (/i = /i2 = 
... = /i„ — > oo), while U amounts to the SWAP oper- 
ation Up (E) crU^ = a ® p. Simple symmetry consider- 
ations show that at the maximum of the initial energy 
tr[i/HO'] the second level is n — 1 fold degenerate, i.e. 
e = £2 = • • • = £«• Denoting 



u — e '^'"'^ oc r2 



we obtain for Qh = Qh(oo) 



Qh{^) = Th\n\ - 



1 



1 + (n- l)u 



(16) 



(17) 



where u is to be found from maximizing the RHS of (17) 
over u, i.e., u is determined via 



1 + (n- l)u + h\u = 0. 



(18) 



Note that in this case W = +oo. In the n 3> 1 
limit we have u = ^^ [1 + o(l)] from (18) and Qh = 



Constrainted maximization 



The case of a finite E in (15) is more complicated. 
We had to resort to numerical recipes of Mathematica 6. 
Denoting {|«H)}fe=i and {\ic)}k=i for the eigenvectors of 
Hu and He, respectively, we see from (11, 12) that W 
and Qh feel U only via the matrix 



a 



ij I kl 



\{^HJcmkHlc)\' 



(19) 
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FIG. 2: The optimal transferred heat Qh versus work W. 
Dashed curves refer to 6* = Tc/Th — 0.9: n — 2 (lower dashed 
curve) and n = 3 (upper dashed curve). Normal curves refer 
to 9 = 0.5: n — 2 (lower normal curve) and n — 3 (upper 
normal curve). 



where u and v in (24, 25) are determined from maximiz- 
ing the RHS of (24) and satisfying constraint (25). 

An important implication of (24, 25) is that Qh{W) is 
an increasing function of W for all allowed values of W: 



Qh{w) > Qh{w') a w> W. 



(26) 



Fig. 2 illustrates this fact. For fixed parameters Tc, 
Th and n, the allowed Ws range from a certain nega- 
tive value — which corresponds to work-extraction from 
a two-temperature system H -I- C — to arbitrary M^ > 0. 
Eq. (26) expresses an intuitively expected, but still non- 
trivial fact that the best transfer of heat takes place upon 
consuming most of the available work. Note that this re- 
sult holds only for properly optimized values of Qh- One 
can find non-optimal set-ups, where (26) is not valid ■^. 



D. Optimization of spontaneous processes 



This matrix is double-stochastic [231 



Y.^.'^^nki-Y.ki'^' 



ij I kl 



= 1. 



Conversely, for any double-stochastic matrix C^- 1 ki there 
is some unitary matrix U with matrix elements Uij i j.; , 
so that Cy I fej = |C/y [feil^ [23]. Thus, when maximizing 
various functions of W and Qc over the unitary U, we 
can directly maximize over the (n^ ~ 1)^ independent 



elements of 



n xn 



^ double stochastic matrix C, 



ij I kl 



This 



fact simplifies the problem. 

Numerical maximization of Qh over all unitaries 
U — alternatively, over all doubly stochastic C matri- 
ces (19) — and energy spectra {f^k}k=2 ^^^ {^fc}fe=2 '^'^^~ 
strained hy W < E produced the following results: 

• The upper energy levels for both systems H and C 
are n — 1 times degenerate [see (3)]: 



fJ-^ fJ-2 



■*"n-) 



£2 



(21) 



• The optimal unitary corresponds to SWAP: 

Up «) crU'' ^<T® p. (22) 

• The work resource is exploited fully, i.e., the max- 
imal Qh is reached for W = E. 

Though we have numerically checked these results for 
n <b only, we trust that they hold for an arbitrary n. 

Denoting by Qh the maximal value of Qh^ and intro- 
ducing from (21) 



V = e-^-^" 



and 



u = e"^"^ 



we have 



Qh 

Th 
W 



In 



{n — l){u — v) 



l + {n - l)v][l + {n - l)uY 
(Inu — 9liiv) (n — l){u — v) 
[l + {n- l)v] [l + {n- l)u] ' 



(23) 

(24) 
(25) 



According to our discussion in section II A, the max- 
imization of transferred heat Qh given by (8) should 



(20) proceed over all unitary operators e 



Id TT 

— s-tl,. 



with iJjnt sat- 
isfying (7) and over the energies {£k}k=2^ {Mfc}JJ=2 ^^ ^ 
and C. This maximization has been carried out along 
the lines described around (20). We obtained that the 
maximal spontaneous heat Qh i^ equal to Qh in (24) 
under condition W — Q: 



sp] 



QhiW^Q). 



(27) 



Thus the optimal spontaneous processes coincide with 
optimal processes with zero consumed work. This re- 
sult is non-trivial, because the class of unitary operators 
with M^ = is larger than the class of unitary operators 
e~"ft ^'"' with Hint satisfying (7). Nevertheless, these two 
classes produce the same maximal heat. 

• Eqs. (26, 27) imply that if the spontaneous heat 
transfer process is already optimal (with respect 
to the junction Hamiltonian) its enhancement 
with help of driven processes does demand work- 
consumption, W > Q. This fact is non-trivial, 
because — as well known from the theory of heat- 
engines — also work-extraction does lead to the heat 
flowing from cold to hot [1, 3]. 

Taking W = Qin (24, 25) and recalhng (23) we get 

p = e, u^v". (28) 



^ The simplest example is a junction, where the free Hamilto- 
nian Hq has a non-degenerate energy spectrum, and thus the 
condition (7) does not hold. There are no proper spontaneous 
processes for this case. Still there can exist a work-exracting 
{W < 0) driven processes that transfer heat from hot to cold. 




w 



FIG. 3: The efficiency x versus work W ior 6 = T^/Tu = 0.5 
and n = 2 (normal curve), n = 10 (dashed curve) and n = 30 
(thick curve). 



The interpretation of (28) is that since there are only 
two independent energy gaps in the system, they have 
to be precisely matched for the spontaneous processes to 
be possible; see in this context the discussion after (7). 

[BP] 



Thus the spontaneous heat Q}^ ' is given as 



2): 



spj 



Tr 



= ln 



(n - l){v'^ - z;o) 



[l + {n-l)vl][l + {n-l)vo] 



(29) 



where vq maximizes the RHS of (29). 



E. How much one can enhance the spontaneous 
process? 

We would like to compare the optimal spontaneous 
heat (29) with the optimal heat Qh{oo) transferred un- 
der consumption of a large amount of work; see (17, 18) 
and recall (26). One notes that for parameters of Fig. 2 
the approximate equality Q;i(oo) k, QhiW) is reached al- 
ready for W/Th < 1. This figure also shows that for the 
temperature ratio 9 = Tc/Th far from 1, the improvement 
of the transferred heat introduced by driving is not sub- 
stantial. It is however rather sizable for 61 ~ 1, because 
here the spontaneous heat (29) is close to zero, while the 
heat Qh{oo) does not depend on the temperature differ- 
ence at all; see Fig. 2 and (17, 18). 



IV. EFFICIENCY 

We saw above that enhancing optimal spontaneous 
processes does require work. Once this is understood, 
we can ask how efficient is this work consumption. The 
efficiency is defined as 



xiw) 



Qh{W) - g^P' 

w 



>o. 



(30) 



where Qh{W) is the optimal heat transferred under con- 
dition that the consumed work is not larger than M^ > 0, 
while Q^ is the optimal spontaneous heat; see (24, 29). 

Note that the two subtracted quantities Qh{W) and Qf^ 
in (30) refer to the same junction H -t- C, but with dif- 
ferent Hamiltonians; see (24, 25). 

For ly -> 0, x{W) increases monotonically and tends 
to a well defined limit x(0); see Fig. 3. 

• For fixed 9 and n, x(0) = x(^ ~^ 0) is the maxi- 
mal possible efficiency at which the enhanced heat 
pump can operate. As seen from Fig. 3, this max- 
imum is reached for 



Qh{W) - Qh{0) -^ +0 and W 



-0, 



(31) 



where we recall that n, Th and Tc are held fixed. 

• There is thus a complementarity between the 
driven contribution in the heat, which according 
to (26) maximizes for W — >■ oo, and the efficiency 
that maximizes under W ^ Q. 

Note from Fig. 4 the following aspect of the maximal 
efficiency x(0): it decreases for a larger n (and a fixed 9). 
This is related to the fact that the optimal spontaneous 
heat Q|^ increases for larger n. 



It is seen from Fig. 3 that 



x{W) < x(0) < 



1 -i 



(32) 



We checked that this upper bound for the efficiency 
(30) holds for all 9 = T^/Th and n. 

It will be seen below that the upper bound y^^ is 
reached in the quasi-equilibrium limit — >■ 1. Note that 
Y3g formally coicides with the Carnot limiting efficiency 
for ordinary refrigerators; see (2). A straightforward im- 
plication of (32) is that enhancing optimal spontaneous 
processes must be inefficient for ^' 0. 

Let us discuss to which extent the bound (32) is similar 
to the Carnot bound (2) for refrigerators. 

0. These two expressions are formally identical. 

1. Recall that (2) is a general upper bound for the ef- 
ficiency of refrigerators that transfer heat against its gra- 
dient. Such a transfer does require work-consumption. 
The same aspect is present in (30), because by its very 
construction the efficiency (30) refers to enhancement 
of the optimal spontaneous process that also demands 
work-consumption. To clarify this point consider a spon- 
taneous process with the transferred heat Qj^^ ■ Let this 
spontaneous process be non-optimal in the sense that 
no full optimization over the Hamiltonians (6, 7) has 



been carried out: Qj^^' < Q^l^ 



[BP] 



ifi ^ ^h ■ This non-optimal pro- 
cess is now enhanced via a work-consuming one. De- 
note by QhiW) > Q^ the transferred heat of this pro- 
cess, where W is the consumed work. Following (30) 



XiO) 




FIG. 4: The maximal efficiency x(0) — xC^ ~ 0) given by 
(33) versus 9 — Tc/T^ for n = 2 (top normal curve), n — 101 
(bottom normal curve), and n = 10"" (dotted curve). Thick 
curve: the efficiency 0/(1 — 6). 



one can define the efficiency of this enhancement as 
X'(W) = [QhiW) - Qt^^]/W. One can now show, see 
Appendix B, that x'(W^) can be arbitrary large for a suf- 
ficiently small (but non-zero) consumed work W. The 
reason for this unboundness is that we consider a non- 
optimal spontaneous process, which can be also enhanced 
by going to another spontaneous process. 

2. We noted above that reaching bound (32) means 
a neglegible enhancement; see (31). The same holds 
for the Carnot bound (2) for refrigerators: operating 
sharply at the Carnot efficiency means that the heat 
transferred during refrigeration is zero; see [19] and ref- 
erences therein. 

3. An obvious point where the bounds (32) and (2) dif- 
fer from each other is that the latter is a straightforward 
implication of the first and second laws of thermodynam- 
ics, while the former is so far obtained in a concrete model 
only. We opine however that its applicability domain is 
larger than this model; some support for this opinion is 
discussed in section V. 



A. Efficiency in the quasi-equilibrium limit ^ ^> 1. 



Let us now find a more explicit expression for x(0) 
using the fact of (31). For a small y^, (25) implies that 
u and V are related as u = v^+x, with x = Ofy^]. Denote 
the RHS of (24) as f[u, v] and the RHS of (25)" as g[u, v 
Expand g[v' 



■ x,v] over x, express x via — and v. 



W 



1 



ngL[v',vr 

and substitute the obtained result into (24): 



Carrying out the optimization of this quantity over v, we 
see that in the first order over the small parameter -y- 
the optimal value reads 

f[vQ,Vo\ + 7^--r-g 7, 

Th gu[v^,vo] 



where vq is defined in (29). The final result for the effi- 



ciency x(0) is I'It, . or 



x(o) = -i 



'ln( — 

.wo 



1 



Vo 



l + {n- l)wg i;o - wg 



.(33) 



Let us now find the asymptotic behaviour of the effi- 
ciency (33) for — s> 1. It is seen from (29) that in this 
limit Vq = a + 0{1 — 9) , where a is found from the tran- 
scendental equation 



In- 



1 _ ^ 1 + (n - l)a 
1 - (n- l)a' 

This implies for x(0) and Q^ ■ 

26 



(34) 



X(0) = — 



-1 



1- (n- l)a 
4Tc (n - 1) a (1 - 6) 



(l-(n-l)a)2 



+ 0(1-0), (35) 
Oiil-ef). (36) 



Thus, in the limit — >■ 1 the efficiency x asymptotically 
reaches the maximal value j^. This reachability is re- 
lated to the fact that the spontaneous heat Q^ becomes 
small; see (36). 

Note that the aymptotics (35) is meaningful only for 



— > 1, where j^ is larger than —1 
Inn 3> 1, (34) implies asymptotically 

a(n- 1) = 1 - 



29 
l-(n-l)a' 



ln[n - 1] 



For 



(37) 



Th g'uiv' 



Thus the limits ^> 1 and n ^ 1 do not commute with 
each other, because the correction term -; — ,^^ ... in (35) 

' 1 — (n — l)a ^ ^ 

will be comparable to the main term j^ . 



V. ENHANCED HEAT TRANSFER IN LINEAR 
NON-EQUILIBRIUM THERMODYNAMICS 

Since the above results were obtained on a concrete 
model, one can naturally question their general validity. 
Here we indicate that these results are recovered from 
the formalism of linear non-equilibrium thermodynamics 
[27-29]. This theory deals with two coupled processes: 
heat transfer between two thermal baths and work done 
by an external field. In contrast to the model studied in 
previous sections, the field is not time-dependent; e.g., 
it can be associated with the chemical potential differ- 
ence [29]. The difference and similarity between set-ups, 
where the work is induced by a time-dependent field, and 



those where the work is done by a constant field is dis- 
cussed in [30] . For mesoscopic models which lead to linear 
non- equilibrium thermodynamics see, e.g., [29, 31]. 

In the linear regime both the driving field and the tem- 
perature difference T^ — Tc are small, so that the heat Qh 
and work W can be linearized [27-29]; see also Appendix 
C. The virtue of this approach is that it is independent 
from the concrete details of the studied system. The role 
of a free parameter — over which the setup can be opti- 
mized — is played here by the phenomenological coupling 
between the heat transfer and work input [27]. 

More specifically, we show in Appendix C that also 
within the formalism of linear non-equilibrium thermo- 
dynamics, enhancing the optimal spontaneous fiow re- 
quires work-consumption. As for the efficiency, we first 
recall the message of (35): for 6 = Tc/Th -^ 1 the max- 
imal efficiency x(0) does not depend on the dimension 
n^ of the junction (provided that the latter is fixed) and 
approaches oc 1/(1 — 9). This asymptotic result is to 
a certain extent universal, because it is also recovered 
within linear non-equilibrium thermodynamics; see Ap- 
pendix C. We stress however that only the asymptotics 
for — > 1 is recovered. Since this is a linear theory re- 
stricted to a small temperature difference and a small 
work input, it is naturally not capable of reproducing 
the full message of the bound (32). For this purpose one 
would probably need a non-linear thermodynamics the- 
ory; see, e.g., [29]. Unfortunately, such theories are not 
so universal (system-independent) as the linear theory. 



VI. SUMMARY 

We started this paper by listing several representative 
examples of enhanced transport and posing two basic 
questions: 

1. When enhancing a spontaneous process does re- 
quire work-consumption? 

2. If such a work-driven enhancement does take place, 
is there a general bound on its efficiency? 



These questions have been answered via a quantum 
model for a heat-transfer junction immersed between 
two thermal baths at different temperatures Tc and Th 
{Tc < Th). The model is defined in section II. We have 
chosen to work with this model of junction, because its 
structure is fiexible enough to allow explicit maximiza- 
tion of the transferred heat over the junction Hamilto- 
nian. We thus can determine the maximal heat trans- 
ferred by the junction. 

Our basic results can be stated as follows. 

1. When the spontaneous heat-transfer is already 
maximized over the junction Hamiltonian, its en- 
hancement does require work-consumption. 

2. The efficiency is defined with respect to the op- 
timal spontaneous heat transfer as the heat incre- 
ment due to ehancement divided over the consumed 
work. This efficiency is shown to be limited from 
above by Tc/{Th — Tc), a bound that is reached for 
Tc —^ Th- For this bound to hold it is essential 
that the efficiency is defined with respect to op- 
timal spontaneous heat transfer. In its turn, the 
very idea of optimality refers to a class of variables 
to be optimized over. For the studied model this 
class is defined by the junctions structure and its 
Hamiltonian. 

The main open problem with these results is whether 
they hold more generally than the studied model. We 
presented a partial evidence — partial because it is re- 
stricted to a small Th — Tc within the linear non- 
equilibrium thermodynamics — that they hold more gen- 
erally. 
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Appendix A: Power of heat transfer 

Recall that the power Q/j/t is defined as the ratio of 
the transferred heat Qh to the cycle duration t. For our 
model T is mainly the duration of the second stage, i.e., 
T is the relaxation time, which depends on the concrete 
physics of the system-bath coupling. 

We discuss two scenarios of relaxation for the present 
model. For the collisional relaxation the target system 
interacts with independent bath particles via successive 
collisions; see [19, 24, 25]. For our purposes the target 
system is H or C that interact with, respectively, the hot 
and cold bath. Each collision lasts a time TcoI, which is 
much smaller than the characteristic time Tbtw between 
two collisions. The interaction Hamiltonian [between the 
target system and a bath particle] is conserved, so that 
there is no energy cost in switching the system-bath in- 
teraction on and off; see [19, 24, 25]. The relaxation time 
T can amount to a few Tbtw^ Tbtw ^ t; see [19]. 

If now the pulse-time d of V{t) is also much smaller 
than Tbtw, one realizes a thermally-isolated pulse (since 
the overlap between the pulse and a collision is negligi- 



ble), while the relaxation time r can be still much smaller 
than any characteristic time of C or H. Thus in this sce- 
nario the cycle time r is finite, and the power of heat 
pumping Qc/t does not vanish due to a large cycle time. 
Within the second scenario the system-bath interac- 
tion is always on, but its magnitude is small [weak- 
coupling]. Now the relaxation time is much larger than 
the internal characteristic time of H and C. Because 
the system-bath interaction is always on, there will be a 
contribution in the work (11) coming from the system- 
bath interaction Hamiltonian. This contribution arises 
even when the conditions for the pulsed regime hold [18]. 
However, within the weak coupling assumption this ad- 
ditional contribution is proportional to the square of the 
system-bath interaction constant and can neglected [18]. 
We stress that this additional contribution does not arise 
within the collisional relaxation scenario, because the 
pulse and collisions are well-separated in time. 



Appendix B 

Let the system H -f C functions as a non-optimal [with 
respect to the Hamiltonian of H + C] spontaneous heat 

pump. The amount of heat transferred per cycle is Qj^ . 
One increases the transferred heat via an external work- 
source acting on H -|- C. Now it is equal to Qh > Qh ■ 
The amount of work consumed per cycle is W . 

The purpose of this Appendix is to show on a concrete 
example that there is no upper bound on the efficiency: 



X 



r 



w 



(Bl) 



It can go to -l-oo for W —^ -1-0. 

Let the systems H and C be three-level systems with 
Hamiltonians [compare with (3)] 

i^H =diag(0,£,£ + ^), i7c =diag(0,^,e + /i) (B2) 

where for future purposes we choose /x > e > 0. The 
initial Hamiltonian (5) reads 

diag[0, e, e -i- ^i, fi, e + ^, e + 2fi, e + fi, 2e + ^, 2e + 2/^] 

The initial state of H + C is given by (2, 3 4): 



P = 



diag[l, a", a^b^] 



1 + a'^ -\-a 
where we denoted 



Bl 



diag[l, b, ab] 
l + b + ab '' 



a = e-^=^ 



b = e-^"''. 



(B3) 



(B4) 



The system H -f C has nine energy levels; three of 
them have equal energy £ + /i. Thus spontaneous pro- 
cesses amount to unitary operations that couple these 
degenerate energy levels to each other, but do not induce 
transitions to non-degenerate energy levels [compare with 



10 



(7)]. Such unitaries will produce double-stochastic ma- 
trices of the following form [compare with (19, 20)] 



/I 
1 



c = 








Cii Ci2 Ci3 

10 

C21 C22 C23 

10.. 

C31 C32 C33 



0\ 














\0 








1 
1/ 



(B5) 



It is not difficult to see that for fixed energy levels (B2) 

the spontanoues process with the largest Q^ amounts to 
Ci3 = C22 ~ C31 = 1, while other c's nullify. This means 
that the spontaneous process interchanges populations of 
the third and seventh energy levels; see (B3). The heat 
transferred per cycle reads 



Ml 



[sp] 



1 



a'^b^ ~ab] (e + fi), 



ah) (1 



2%") 



(B6) 

(B7) 



Now the external field acts on the system H + C en- 
hancing the heat transfer, i.e., increasing Q^ ■ We pos- 
tulate that this action amounts to SWAP transformation. 
We obtain for the heat, work and efficiency (Bl) [com- 
pare with (12, 11)]: 



1 



{ey 



-6H 
ab 



ba\a-b'^)] 
ba'^ib^ -a)]}, 



0] 



+(e + ^l) [(ab) 

W^-ifi-e)[a'' -b- ba^b" - a)] 

It is seen that when e — > /i [which means a 

the difference Qh — Qh ^^ positive and finite, while the 
work W is positive, but turns to zero as W cc fx — e. 
Thus if one defines efficiency as x' = (Qh — Qh )/W, 
it will turn to infinity for £ — !> /x. The divergence of x' 
is ultimately due to the fact that the spontaneous heat 
was maximized over Hamiltonian (6) only partially, i.e., 
it was not maximized over the energy levels of H and C. 



work done by an external time-independent force. Thus 

Xi^Pc-Ph = ^^^Y^^ ^2 = ^, r = V7)j;,(C2) 

where / is the force responsible for the work [28]. Once 
we are in the linear regime over small parameters Th — Tc 
and /, one can substitute T in (C2) by Tc or Th] the 
choice of T is conventional. 

In essence, linearity means that the state which sup- 
ports the currents is not far from equilibrium [28]. The 
basic postulate of this formalism is linear relations be- 
tween currents and forces [28]: 



Ji — LiiXi + L12X2 



J2 — L21X1 + _L22^2- (C3) 



where the kinetic coeffcients Lik do not depend on Xi due 
to assumed linearity of the overall process. The kinetic 
coeffcients L12 and L21 quantify the coupling between 
the two processes. 

The statement of the second law relevant for this com- 
posite linear process amounts to the positivity of entropy 
production diS/dt [28] 



d^ 
dt 



Xi J I -\- X2 J 2 



2 

E 

i,fe=l 



UkX.Xk > 0. (C4) 



As a consequence of the time-invariance of the underlying 
microscopic theories (i.e., classical or quantum mechan- 
ics) , the matrix of kinetic coefficients is symmetric [28] 



^21- 



(C5) 



We now assume that some work is dissipated, J2 > 0, 
for enhancing the heat flow. The spontaneous heat trans- 
fer corresponds to no coupling between the processes: 
L12 — L21 — 0. 

In calculating the efficiency of enhancement we shall 
follow the same strategy as in sections III and IV: first 
we shall maximize the transferred heat Ji under a fixed 
amount of work W = / J2 = TX2 J2. Then the efficiency 
will be defined as in (30). The maximization variables 
are /, Ln, L22 and L12 = -/j21- The temperatures Tc and 
Th are held fixed; see also (C2). Write Ji as 



Appendix C: Heat-transfer enhancement from linear 
non-equilibrium thermodynamics 

The formalism of linear non-equilibrium thermody- 
namics starts by introducing currents of physical quanti- 
ties Ji (e.g., currents of energy or mass) and the respec- 
tive conjugate forces Xi [27, 28]. We introduce two such 
currents and forces: 



^1 , J2, A"i , X2 



(CI) 



Ji will refer to heat flowing from the high-temperature 
bath at temperature Th = i/l3h to the lower one at tem- 
perature Tc = l//3c, while the second current J2 is the 



Ji — LiiXi 



W X| 

TXi Xi 



(C6) 



During the maximization we should keep Ln confined 
by some upper limit Cn; otherwise Ji will not be finite. 
Eq. (C5), which should hold for arbitrary Xi and X2^ 
implies ^22 > 0. Recaling that also Xi > 0, we see that 
Ji is maximized for Ln — £11 and L22 ~ 0: 



Ji — jCiiXi 



W 
7% 



(C7) 



Thus for enhancing the optimal spontaneous heat CnXi 
we need W > (work-consumption). Subtracting from 



11 

J\ the spontaneous contribution C\\X\ and dividing over which for a smah temperature difference T^ — T^. (recall 
the consumed work W , we get for the efficiency: that this is the applicability domain of the considered 



1 T 



linear theory) coincides with (35). 



